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SUMMARY 

Pitot rakes are often used to measure mass flow rate In circular or 
annular ducts. Often the rakes are area weighted and a simple summation Is 
used to determine the average velocity. Errors In flow rate measurement are 
Inherent In this technique because of the discretization of the velocity pro- 
file. The error decreases as the number of tubes on the rakes Increases, and 
resolution of the velocity profile Improves. A study was conducted to deter- 
mine the error In measuring mass flow rate with pitot rakes In an annulus. 

The Ideal flow rate was determined by using a unique semlemplrlcal analysis 
for fully developed, turbulent flow. The velocity profile obtained from this 
analysis was Imposed on the pitot rake, and an area-weighted summation was 
used to determine the flow rate that the rake would Indicate. Results In terms 
of flow coefficient, or the ratio of Ideal to Indicated flow rate, ranged from 
0.903 for one probe placed at a radius dividing two equal areas to 0.984 for a 
10-probe area-weighted rake. Flow coefficients were not a strong function of 
annulus hub-to-tlp radius ratio for rakes with three or more probes. 


INTRODUCTION 

The measurement of airflow rate Is of primary Interest In the testing and 
development of aircraft propulsion systems. Several methods of measuring flow 
rate are available Including choked-exlt devices, venturis, and bellmouths. 
These methods generally require the Installation of additional hardware for 
the specific purpose of measuring flow rate. Where Installation of such hard- 
ware Is Impractical, other means of measuring flow rate must be devised. For 
Inlet testing there Is usually an array of total pressure rakes at the diffuser 
exit station to measure total pressure recovery and distortion. Such rakes 
are generally area weighted and may also be used to measure flow rate. Each 
ring (probe In each rake at a common radius) In the array Is assigned an area 
of the circular or annular duct. If the rakes are area weighted, these are 
equal areas. This area combined with the measured ring total pressure, a 
representative static pressure, and the total temperature can be used to cal- 
culate a mass flow rate for that ring. Subsequent summation of all ring flow 
rates results In an estimate of the total mass flow rate. Because of the dis- 
crete nature of the measurements significant error exists In the Indicated 
value of flow rate, especially In high velocity gradient regions such as boun- 
dary layers. Obviously the larger the number of probes on the rakes, the finer 
the resolution of the velocity profile and the smaller the error. To account 
for this discrepancy between the measured and Ideal flow rates, experimentally 
determined flow coefficients applicable to a particular rake geometry and 
annulus hub-to-tlp radius ratio have normally been used (ref. 1). However, 





since the flow coefficient varies with the number of probes on the rakes and 
the hub-to-tlp radius ratio of the annulus, a semlemplrlcal procedure for 
determining an applicable flow coefficient would be desirable in many cases. 

This Investigation was conducted to determine the error In the rake mass 
flow rate measurement for turbulent flow In an annulus due to the discretiza- 
tion of the velocity profile and the variation of this error with the number 
of probes on an area-weighted rake. Also Investigated was the effect of annu- 
lus hub-to-tlp radius ratio on the error. 

A semlemplrlcal method for determining the velocity profile In fully 
developed, turbulent flow In an annulus Is presented. The Integral of this 
velocity profile Is compared with summation of velocities as measured by an 
area-weighted rake subjected to the same velocity profile. Results are pre- 
sented as flow coefficients, or the ratio of the Integrated profile to the 
summed, for area-weighted rakes of up to 10 probes and a range of annulus hub- 
to-tlp radius ratios. 


SYMBOLS 

A^ area assigned to probe 1, m^ 

A n Integration constant for region n 
B n constant for region n 

Cp flow coefficient 

C n constant for region n 

D n constant for region n 

E n constant for region n 

F(n) function of region n 

F(x) function of variable x 

I number of probes on rake 

K n dimensionless distance from origin to Inner radius of region n, r/R^ 

Kp dimensionless distance from hub surface to Inner radius of region n, 

(r - R h )/(Rt - Rh) 

m n ratio of turbulent to laminar viscosity In region n 
N number of regions In semlemplrlcal method 

p static pressure gradient In flow direction, Pa/m 

r distance from origin In radial direction, m 
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Rh 

Rt 


annulus hub radius, m 

annulus tip radius, m 

u dimensionless distance from origin In radial direction, r/Rt 

v velocity In z-dlrectlon, m/sec 

v max maximum velocity In z-dlrectlon, m/sec 

V n velocity In z-dlrectlon at Interface of regions n and n - 1, m/sec 

^(8-) coefficient of laminar viscosity, Pa sec/m 1 2 

v (t) coefficient of turbulent viscosity. Pa sec/m 2 

(9.1 2 

rz laminar transport of z-momentum In r-dlrectlon. Pa/m 

(tl 2 

t' turbulent transport of z-momentum In r-dlrectlon, Pa/m 

Subscripts: 
h hub surface 

1 Individual probe on rake 

N outermost region In semlemplrlcal method 

n Individual region In semlemplrlcal method 

r radial direction 

t tip surface 

z axial direction 

Superscripts: 

(l) laminar 

(t) turbulent 

APPROACH 

Derivation of Semlemplrlcal Velocity Profiles 

To develop the turbulent velocity profiles, a semlemplrlcal approach Is 
used under the following assumptions: 

(1) Two-dimensional (axlsymmetrlc) flow 

(2) Steady, fully developed flow 
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(3) Incompressible flow 

(4) Linear pressure gradient In flow direction 

Figure 1 depicts the flow situation and the coordinate system. The time- 
smoothed z-component of the momentum equation under the preceding assumptions 
becomes 



Before equation (1) can be Integrated to yield velocity profiles, an 
additional equation Is needed to relate the laminar and turbulent flux of 
z-momentum In the radial direction to the radial velocity gradient. To develop 
this relation, the annular cross section Is divided Into an arbitrary number 
of annular regions N. Within each region It Is assumed that the negative of 
the radial velocity gradient Is proportional to the sum of the laminar and 
turbulent momentum fluxes. Further the constant of proportionality In a 
particular region n Is the sum of the laminar viscosity and a turbulent 
viscosity that varies from region to region. In equation form 


.00 + t (t) 


rz 


rz 


-(• 


») , u <t)\av 

v n /dr 


( 2 ) 


Expressed In terms of the ratio of turbulent to laminar viscosity, which Is 
constant In a particular region n, equation (2) becomes 


(a) . (t) 


(it), 


dv 


where 


rz + T rz v ^ + m n^ dr 

..(*) 

n 


m " = Jit) 


Combining equations (1) and (3), there results for a region n 

P = 

Integrating over the radial direction yields after rearrangement 

1 


1 L L (*•)/! , m s dv 

F dF L (1 dF 


dv 

dr 


v U) 0 * m n ) 


' A \ 

pr _n 

,2 r 


(3) 

(4) 

(5) 

( 6 ) 


where A n Is the Integration constant for region n. Defining K n as the ratio 
of the Inner boundary radius of region n to the annulus tip radius Rt and 
defining V n as the velocity at the Inner boundary radius of region n. 
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equation (6) can be Integrated from the Inner boundary of region n to an 
arbitrary radius within the region: 



Finally for any region n 


v(r) = V r 


1 


y (8,) (l ♦ m n ) 


pR t r 




- K n] * A n >" (VrJ 


( 7 ) 


(8) 


Note that K-| Is equal to the hub-to-tlp radius ratio and Kn+t Is equal to 1. 
At the Interface of two adjacent regions the solutions are "spliced" together 
by equating the velocity gradients. Writing equation (6) for regions n and 
n + 1 at a common radius K n+ -j and equating the two resulting expressions 
yields 


1 + m_ 


[PK n R + 
K n+lt 


K n+l R t, 


1 


'P K n+l R t 


n+1 


1 + m, 


n+1 


K n+l R t 


(9) 


Solving for results In a recursion relation 


n+1 


"l + m , 
n+1 | 

^♦l R t fl n ^ 


1 + m 
L n 

V 2 K n+1 R tJ 

1 

CVI 

1 


K n+l R t 


( 10 ) 


Velocities at the Interfaces between regions can be computed by evaluating 
equation (8) at r = K n+ -|Rfc: 



Applying the no-slip boundary condition at the tip radius R, equation (11) 
becomes 


1 


v (!l) (l + m N ) 


PR 




2 

N+1 



( 12 ) 
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Now V n can be successively replaced by using equation (11) until an equation 
of the form 



results. Since Is the velocity at the Inner boundary of region 1, It Is 
the velocity at the hub surface and Is zero. Hence 

N 

O VE F(n) (14) 

n=l 


Equation (14) along with the N - 1 equations obtained from the recursion 
relation for A n (eq. (10)) can be arranged Into the following form: 


C 1 A 1 + C 2 A 2 + * * • + C N-1 A N-1 + C N A N = B 
E, A, - A„ = -D, 


VI "2 

E 2 A 2 - A 2 


= -D, 


(14) 

( 10 ) 

(10) 


where 


E N-lVl - a n ■ - d n-i <10 > 


B P R t Y - < 

" 4 4v* ra " 


c . -4- 

" 1 * m n \ K nV 


P(K nt 1 R t ) ] * Vi 
" ' 2 \ ’ * \ 


1 + m 


n+1 


1 + m 


n 


(15) 


(16) 


(17) 

(18) 
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Cramer's rule can be used to solve for the Integration constant In region 1 
(A-]), after which repeated application of equation (10) yields all other A n 's. 
The velocity profile for the entire annulus Is now computed by using equation 
(8), beginning at the hub surface (where V-| = 0) and progressing to the 
outer radius of region 1. At the outer radius of region 1 a solution for 
V 2 results, and equation (8) Is applied again for region 2. In a similar 
manner the solution progresses across the annulus to the tip surface, where 
the no-slip boundary condition has previously been satisfied In equation (14). 

Values of turbulent-to-lamlnar viscosity ratio m n and region bound- 
aries, as well as the number of regions used N In this Investigation are 
depicted In figure 2. These constants were found to best fit available experi- 
mental velocity profile data over a wide range of hub-to-tlp radius ratios 
with only the region 5 viscosity ratio ms varying with hub-to-tlp radius 
ratio. In figure 3 nondlmenslonal velocity profiles obtained by this procedure 
show good agreement with the experimental data of reference 2. 


Calculation of Flow Coefficients 

For the purpose of this analysis flow coefficient Is defined as the Ideal 
flow rate divided by the measured or Indicated flow rate. The Ideal flow rate 
Is determined by using the semlemplrlcal velocity profile. The Indicated flow 
rate Is determined by using an area-weighted summation of velocity values taken 
from the semlemplrlcal profile at the rake probe locations. Under the assump- 
tion of Incompressible flow the flow coefficient reduces to the following: 

1 

2 / v(u)u du 


<’ - K ?> E Vi 

1=1 

where v^ Is the velocity at probe location 1, A^ Is the annular area 
fraction associated with that probe, and I Is the number of probes on the 
rake. If the rake Is area weighted, all A^ are equal. Solutions to 
equation (19) were obtained by using the Fortran computer program listed In 
the appendix. The program evaluates the numerator by applying a sixth-order 
numerical Integration technique known as Weddle's method to equation (8). 

Inputs to the program are annulus hub-to-tlp ratio and the number of probes on 
the rake. A subroutine automatically locates the probes at area-weighted radii 
although provisions exist to Input other rake geometries through an Input data 
set. Turbulent-to-lamlnar viscosity ratios, region boundaries, and the number 
of regions used In conjunction with equation (8) are Input through a separate 
Input data set. All results presented were obtained by using the values In 
figure 2. 




09) 
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RESULTS 


Figure 4 depicts the differences between the Ideal velocity profile and 
the discretized profile that Is summed In obtaining the Indicated flow rate. 

For all probes compensating errors occur. However, probes near the hub and 
tip surfaces will clearly overpredict the flow rate because of the no-slip 
condition at these surfaces. The discretized profiles more closely approximate 
the Ideal profile as the number of probes Increases. Figure 5 presents the 
flow coefficients obtained with the Fortran program for area-weighted rakes 
with one to 10 probes over a range of annulus hub-to-tlp radius ratios. 


CONCLUSIONS 

An Investigation was conducted to determine the error In the rake mass 
flow rate measurement for turbulent flow In an annulus due to the discretiza- 
tion of the velocity profile and the variation of this error with the number 
of probes on an area-weighted rake. The following conclusions were drawn: 

1. The semlemplrlcal method presented for determining fully developed, 
turbulent velocity profiles In an annulus agreed adequately with experimental 
data. 


2. Flow coefficients ranged from 0.903 for one probe placed at a radius 
dividing two equal areas to 0.984 for a 10-probe area-weighted rake. 

3. Flow coefficients were not a strong function of annulus hub-to-tlp 
radius ratio for rakes having three or more probes. 
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APPENDIX - COMPUTER PROGRAM LISTING 


OOOOlOO 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

0001200 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 

0002300 

0002400 

0002500 

0002600 

0002700 

0002800 

0002900 

0003000 

0003100 

0003200 

0003300 

0003400 

0003500 

0003600 

0003700 

0003800 

0003900 

0004000 

0004100 

0004200 

0004300 

0004400 

0004500 

0004600 

0004700 

0004S00 

0004900 

0005000 

0005100 

0005200 

0005300 

0005400 


CXXXXXDONUT. SOURCE. TAPVA 

c 

C THIS PROGRAM COMPUTES SEMI-EMPIRICAL VELOCITY PROFILES FOR FULLY-DEVELOPED TURBULENT 
C FLOW IN AN ANNULUS. THESE PROFILES ARE COMPARED TO DISCRETIZED PROFILES AS INDICATED 

C BY A PITOT RAKE SUBJECTED TO THE SEMI-EMPIRICAL FLOW FIELD. FLOW COEFFICIENTS ARE 

C CALCULATED BY DIVIDING THE 'ACTUAL' FLOWRATE BY THE INDICATED FLOWRATE. THE 'ACTUAL* 

C FLOWRATE IS OBTAINED BY A NUMERICAL INTEGRATION OF THE SEMI-EMPIRICAL PROFILE. THE 

C INDICATED FLOWRATE IS OBTAINED BY AN AREA-WEIGHTED SUM OF VELOCITIES OBTAINED FROM 

C THE SEMI-EMPIRICAL VELOCITY PROFILE AT THE RAKE PROBE LOCATIONS. RESULTS ARE OUTPUT 

C GRAPHICALLY. 

C 


DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 


CAPPA(10),VISCO(10) 

VISCOL.PRESS.RADIUS 

B(10),C<10),D1,D2,D3,DC10),E(10),BSUM,ATERM.ANUM,ADEN,AC10) 

VI2000), VINTER (10), RATIO, VREF,H(10),HSUB(10), VS UM, SUBS UMC10), PI, P2,P3,RV,C0EFF(7) 


DIMENSION NINT(10),VEL(2000).RAD(2000) 

DIMENSION AREA( 100 ),TLOC( 100), DISCI 100 ),PLOC( 100 ) ,TVEL( 100) ,PVEL(100) 
DIMENSION PFL0UC2) 

DATA COEFF/1.,5.,1.,6.,1.,5.,1./ 


DIMENSION XVARS(lO).YVARSaO) 

DIMENSION IVARS(IO) 

C 

C 

CXXXXXINPUT REGION GEOMETRY AND FLOW VARIABLES 
C 

READC5.1000) NREG 
READ( 5 , 1500 ) ALPHA 
DO 100 1=1, NREG 

READ(5, 1100) CAPPA(I),NINT(I),VI5CO(I) 
CAPPAC I )=ALPHA + ( 1-ALPHA > XCAPPA( I ) 

100 CONTINUE 
C 

READC5, 1200) VISCOL 
READC 5 , 120 0 ) PRESS 
READ! 5,1200) RADIUS 
C 

CAPPA(NREG+1)=1. 

C 

C 

CXXXXXCOMPUTE INTEGRATION CONSTANT FOR REGION ONE 
C 

DO 200 N=1,NREG 
C 


B(N)=PRESSXRADIUSXX2X(CAPPA(N+l)XX2-CAPPA(N)XX2)/(4X(l+VISCO(N)>) 

C(N) = CDLOG(CAPPA(N+1)/CAPPA(N)))/(1+VISCO(N)) 

Dl=(l+VISCO(N+l))/(l+VISCO(N)) 

D2=PRESSX(CAPPACN + 1)XX2)X(RADIUSXX2V2 
D(N) =D2X( Dl-1 ) 

E(N)=D1 


200 CONTINUE 


0005500 

0005600 

0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006400 

0006500 

0006600 

0006700 

0006800 

0006900 

0007000 

0007100 

0007200 

0007300 

0007400 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008400 

0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 


C 

C 

CXXXXXSOLVE NXN MATRIX FOR AC1) AND COMPUTE REMAINING A(N)'S 
C 

BSUM=0 

DO 220 1=1, NREG 
BSUM=BSUM-B ( I ) 

220 CONTINUE 
C 

ATERM=D( 1 ) 

ANUM=BSUM 
DO 240 1=2, NREG 
AHUM=ANUM-C( I )XATERM 
ATERM=ATERMXE(I)+D(I) 

240 CONTINUE 
C 

ATERM=1 

ADEN=0 

DO 260 1=1, NREG 

ADEN=ADEN+C(I)XATERM 

ATERM=ATERMXE(I) 

260 CONTINUE 
C 

AC1)=ANUM/ADEN 

C 

DO 300 N=2,NREG 
ACN)=D(N-1)+ECN-1)XACN-1) 

300 CONTINUE 
C 
C 

CXXXXXNUMERICALLY INTEGRATE VELOCITY PROFILE FROM CAPPAC1) TO 1 

C 

C 

V5UM=0 . 

VINTERC 1 ) =0 . 

JSTART = 1 
C 

DO 600 N=1 , NREG 


C 
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0009400 SUBSUflt N) =0 . 

0009500 H(N)=CCAPPACH+1)-CAPPA<N))/NINT<N> 

0009600 HSUB<N > =H( N)/6 . 

0009700 RATIO=CAPPA(N) 

0009800 C 

0009900 NSTOP=NINT<N> 

0010000 DO 500 NSUB=1,MST0P 

0010100 C 

0010200 JST0P=JSTART*6 

0010300 K=1 

0010400 C 

0010500 DO 400 J=JSTART, JSTOP 

0010600 C 

0010700 Pl=l ./( VISC0LM(1+VISC0(N) ) ) 

0010800 P2=PRESSMRADIUSmM2m(RATI0MM2-CAPPA<N)MM2)/4 

0010900 P3=A(N)MDL0G(RATI0/CAPPA(N>) 

0011000 V(J)=VINTER<H>-P1M<P2+P3) 

0011100 RV=V( JJMRATIO 

0011200 C 

0011300 RAD(J)sRATIO 

0011400 SUBSUM(N)=SUBSUM(N)+COEFF<K)MRV 

0011500 K=K+1 

0011600 RATIO=RATIO+HSUB(N) 

0011700 C 

0011800 400 CONTINUE 

0011900 C 

0012000 RATIO=CAPPACN)+NSUBMH<N) 

0012100 JSTART=J 

0012200 C 

0012300 500 CONTINUE 

0012400 C 

0012500 SUBSUM(N)=.3NHSUB(N)NSUBSUM(N) 

0012600 VSUM=VSUM+SUBSUM(N> 

0012700 VINTER(HU)=V(J) 

0012800 C 

0012900 600 CONTINUE 

0013000 C 
0013100 C 

0013200 CM M MM M COMPUTE POINT OF MAXIMUM VELOCITY IVREF), AND VAVO 
0013300 C 

0013400 N=(NREGFl)/2 

0013500 C 

0013600 RATI0=DS<JRT(-2mA<N)/(PRESSMRADIUSMM2)) 

0013700 C 

0013800 Pl = l./<VISCOLMd+VISCO<N)>> 

0013900 P2=PRESSMRADIUSMM2M(RATI0MM2-CAPPA(N)MM2)/4 

0014000 P3=ACN)MDLOG(RATIO/CAPPA(N>> 

0014100 VREF=VINTER(N)-P1M(P2*P3) 

0014200 C 

0014300 VAVG = 2MVSUM/( VREFM(l-CAPPAt 1JMM2) } 

0014400 C 
0014500 C 

0014600 CMMMMMPLOT THEORETICAL PROFILE USING POIHTS OF INTEGRATION 
0014700 C 

0014800 DO 650 1=1, JSTOP 

0014900 VEL(I)=V(I)/VREF 

0015000 650 CONTINUE 

0015100 C 

0015200 XVARS<1)=9 

0015300 XVARS ( 2 ) =6 , 

0015400 XVARS(3)=0. 

0015500 XVARS(4)=0. 

0015600 XVARS 1 5) =1 . 

0015700 XVARS(6)=5. 

0015800 XVARS(7)=1. 

0015900 XVARS (8) = 6 . 

0016000 XVARS(9)=CAPPA(1) 

0016100 CALL XAXISCl.,1.. XVARS) 

0016200 C 


0016300 YVAR5(1)=9 

0016400 YVARS(2)=8. 

0016500 YVARS(3)=90. 

0016600 YVARS(4)=CAPPA(1) 

0016700 YVARS(5)=1. 

0016800 YVARS ( 61=5. 

0016900 YVARS ( 7 ) = 1 • 

0017000 YVARS(8)=2. 

0017100 YVARS(9)=0. 

0017200 CALL YAXISU 1 . .YVARS) 

0017300 C 

0017400 CALL CORNER(l) 

0017500 C 

0017600 IVARS(1)=2 

0017700 IVARS(2)=JSTART 

0017800 C 

0017900 CALL GPLOT(VEL.RAD.IVARS) 


0018000 C 

0018100 CALL AVRAD(VEL.RAD.JSTART,VAVG,RADH,RADT> 

0018200 C 

0018300 CMMMMMANALYZE RAKE PROFILE 

0018400 CMMMMMINPUTS ARE NUMBER OF TUBES, RADIAL LOCATIONS 
0018500 CMMMMMAND ANNULAR AREA FRACTION ASSUMED FOR EACH TUBE 
0018600 C 

0018700 TSUM = 0 , 

0018800 READ(6,1300) NTUBE 


0018900 C 
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00X9000 

0019100 

0019200 

0019300 

0019900 

0019500 

0019600 

0019700 

0019800 

0019900 

0020000 

0020100 

0020200 

0020300 

0020900 

0020500 

0020600 

0020700 

0020800 

0020900 

0021000 

0021100 

0021200 

0021300 

0021900 

0021500 

0021600 

0021700 

0021800 

0021900 

0022000 

0022100 

0022200 

0022300 

0022900 

0022500 

0022600 

0022700 

0022800 

0022900 

0023000 

0023100 

0023200 

0023300 

0023900 

0023500 

0023600 

0023700 

0023800 

0023900 

0029000 

0029100 

0029200 

0029300 

0029900 

0029500 

0029600 

0029700 

0029800 

0029900 

0025000 

0025100 

0025200 

0025300 

0025900 

0025500 

0025600 

0025700 

0025800 

0025900 

0026000 

0026100 

0026200 

0026300 

0026900 

0026500 

0026600 

0026700 

0026800 

0026900 

0027000 

0027100 

0027200 

0027300 

0027900 

0027500 

0027600 

0027700 

0027800 

0027900 

0028000 

0028100 

0028200 

0028300 

0028900 

0028500 

0028600 

0028700 


DO 700 I=1,NTUBE 
READ(6,1900) TLOC( I ) , AREA! I) 

700 CONTINUE 
C 
C 

CXXXXXCOMPUTE points of discontinuity 

c 

AT0T=1-CAPPA(1)XX2 

DISC(1) = CAPPA(1) 

c 

ISTOP=NTUBEH 
DO 750 1=2, I5TOP 

DISCI I ) -SORT (CAREA(I-1)xAT0T)+DISC(I“1)xx2) 

750 CONTINUE 
C 

c 

c 

CXXXXXCOMPUTE VELOCITY AT EACH TUBE LOCATION AND DO WEIGHTED SUM 
C 

DO 800 J-1,NTUBE 
C 

DO 850 NM.NREG 

IF(CAPPA(NU).GE.TLOC!J))GO TO 875 
850 CONTINUE 
875 CONTINUE 
C 

RATIO=TLOC(J) 

C 

P1=1./(VISCOLX(1+VISCO(N))) 

P2=PRESSXRADIUSXX2x(RATIOXX2-CAPPA(N)xx2)/9 

P3=A(N)XDLOG(RATIO/CAPPA(H)> 

V( J) = < VINTER! N) -PI x(P2 + P3 > )/VREF 
TSUM=T5UM+AREA(J)XV(J) 

C 

800 CONTINUE 
C 
C 

CXXXXCOMPUTE FLOW COEFFICIENT 
C 

CFLOU=VAVG/TSUM 

C 

C 

CxxxxxpLOT 'ASSUMED* PROFILE, WITH TUBE LOCATIONS 
C 

NP0INT=HTUBEX2 

C 

DO 900 I=1,NTUBE 
II=IX2 

PVEL ( I I )=V( I ) 

PVEL(II“1)=V(I) 

900 CONTINUE 
C 

PLOC(l)=CAPPA(l) 

DO 950 I=1,NTUBE 
II=IX2 

PLOC(II)=OISC(I+l) 

PLOC(II+l)=DISC(I+l) 

950 CONTINUE 
C 

IVARS(2)=NP0INT 

C 

CALL GPLOTCPVEL ,PLOC, IVARS) 

C 

C 

DO 975 I=1,NTUBE 
TVEL! I ) =V( I ) 

975 CONTINUE 
C 

IVARS(1)=6 
IVARS(2)=NTUBE 
IVARS(3)-3 
IVARS(9)=62 
IVARS ( 5 ) = 1 
IVARSC61-20 
C 

CALL GPLOTtTVEL.TLOC, IVARS) 

C 

C 

CALL NUMBER(9.CFL0W.6.9,PFL0W) 

CALL CHARSC6,PFLOW,0,9.8, .5,15) 

CALL CHARS! 12, 'FLOW COEFF = ' , 0 , 3 . , .5, 15 ) 

C 

PRIHT1500.CFLOW 
CALL DISPLA! 1 ) 

C 

STOP 

C 

C 

CXXXXXFORMATS 

c 

1000 F0RMATII2) 

1100 FORMAT!F6.9,2X,I2,2X,F10.5) 

1200 FORMAT (Dll. 9) 

1300 FORMAT (13) 

1900 FORMAT(F10.5,2X,F10.5) 

1500 FORMAT ( F10 .5) 

C 

END 
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